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THE STRAIGHT LINE SOLUTIONS OF THE 
PROBLEM OF If BODIES* 

By F. R. Moulton 

1. Introduction. In an elegant prize memoir,! Lagrange proved that 
there are, for any three finite masses attracting one another according to the 
Newtonian law, four distinct configurations such that, under proper initial pro- 
jections, the ratios of the mutual distances remain constants. In all the con- 
figurations the bodies describe similar conic sections with respect to the center 
of mass of the system, the simplest case being that in which the orbits are cir- 
cular. In three of the four solutions the masses lie always in a straight line, 
and in the fourth they remain at the vertices of an equilateral triangle. 

In this paper two closely related problems are treated : 

I. The number of straight line solutions is found for n arbitrary positive 
massesj ; that is, the ratios of the distances ai-e determined so that under 
proper initial projections the bodies will always remain collinear. This is 
the direct generalization to the problem of n bodies of Lagrange's straight line 
solutions. The method of Lagrange is not suitable for treating the general 
case. Below the problem is formulated so that it is solved by finding the total 
number of real solutions of n simultaneous fractional algebraic equations. 
From the mathematical standpoint the interest in this part of the paper centers 
in the algebraical problem of finding the number of real solutions of a formid- 
able looking set of simultaneous equations. 

II. Tlie second problem is that of determining, when possible, n masses 
such that if they are placed at n arbitrary collinear points, they will, under proper 
initial projection, always remain in a straight line. If n is even and the linear 
dimensions of the orbits are given, it is proved that the n masses are in gen- 
eral uniquely determined ; and that if n is odd the coordinates of the n points 
must satisfy one algebraic relation, after which, choosing any one of the 

* The greater part of this paper was written in 1900 and was presented to the Chicago 
Section of the Am. Math. Soc, Dec. 28, 1900. See Brill. Am. Math. Soc, vol. 7 (1900-1901), 
p. 249. It has been recast and put in form for publication by the author as a Research 
Associate of the Carnegie Institution of Washington. 

t Lagrange, Collected Works, vol. 6, pp. 229-324. 

t Treated briefly by Lehmann-Filhfes in Astronomische Nachrichten, vol. 127(1891), 
no. 033. 
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masses arbitrarily, the remaining n— 1 are uniquely determined, 
problem is in a certain sense the converse of the preceding. 



This 



I. Determination of the Positions when the Masses are given. 

2. The Equations Defining the Solutions. Let the origin of 
coordinates be taken at the center of gravity of the system. This point and 
the line of initial projection of any mass determine a plane. All the other masses 
must be projected in this plane, for otherwise they would not be collinear at the 
end of the first element of time. All the bodies being initially in a line and pro- 
jected in the same plane, they will always remain in this plane. Conse- 
quently if solutions exist in which the n masses are always in a straight 
line, the orbits are plane curves. 

Let the plane of motion be the f »/ — plane. Let the masses be denoted by 
m^, wij, • • •, »ra,„ and their respective coordinates by (^i, t/j), (^^, rj^), . . ., 
(?«> Vn)' Then, choosing the units so that the Gaussian constant is unity, 
the differential equations of motion are 
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Equations (1) admit the integral of areas 
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In case the n bodies remain collinear the line of the resultant accelera- 
tion to which each one is subject always passes through the origin. Therefoi'e 
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in collinear solutions we have, by the law of areas, for each body separately 

miii -^ =Ci , (i = l, . . -, n). 
But when the bodies remain collinear we have also 



ddi 


de^ 
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dt ~ 


dt 


~ dt' 


from which it follows that 






(3) ? = 


Imjc, 


= an. 



where the a^ are constants. That is, if any collinear solutions exist the ratios 
of the distances of the bodies from the origin are constants, and it easily 
follows from this that the ratios of their mutual distances arc also constants. 
They are therefore of the Lagrangian type. 

If the n masses remain collinear, the ratios of their coordinates are con- 
stants, being equal to the ratios of their distances from the origin. There- 
fore in all collinear solutions we have 

(^) h = Xi^, Vi = ^iV, (i = 1, • • •, n), 

where the Xi are constants. Substituting in equations (1) we have, as neces- 
sary conditions for the existence of the collinear solutions, 

n 

d^~~ l^i x,[{x,-xjr-r^ -^ , J ^ t, C* _ 1, . . ., n), 
(5) { d^T, ^ mj(Xi - xj) V 

di^-- fy Xi[{Xi-xjyy/^ r*' 



r = y/f^ + iK 

In order that f and ij as^defined by their initial^ values and equations (5) 
shall be the same for all values of i, the coefficients of !/>■* and i]/r^ must be 
set equal to a constant independent of i. Letting — w'^ represent this constant 
and Vij = ^ {Xi — Xj)'^, these conditions, which are sufficient as well as neces- 
sary for the existence of the collinear solutions, become 
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m^jxi - X.2) mgjxi - X3) m„{xi - x„) 

U + -3 + -3 +••• + -5 =coXi, 

'12 '13 '1» 

[ml(x^ -Xi) ml{x2 - xl) tn„(x^ - x„) 



(6) 



OT,(x„-a;i) mi(x„-Xi) ms(x^ - Xg) 

^s 1 3 r Ts !-••• + 



»il <i 'n3 



= <<)'X„ 



It will be shown that aP must be positive in order that real solutions of 
these equations shall exist. Suppose the notation is chosen so that in any 
solution Xl < • ■ ■ < x„- Then the terms of the left member of the last equa- 
tion are all positive. Therefore «' > unless a;„ < 0, which is impossible 
since the center of gravity relation mj |i + • • • 4- m„ f„= holds true. 

For every set of real values of Xi, • • • ,x„ satisfying equations (6) the 
solutions of (5) are the same for all i, and these solutions substituted in (4) 
give the coordinates in the collinear configurations. 

Since equations (5) have the same form as the differential equations in 
the two-body problem, it follows that in the collinear solutions the orbits are 
similar conic sections. In case the orbits are ellipses the coefficient of 
— l/r-^ and — 7//r^ is the product of the cube of the major semi-axis of the orbit 
and the square of the mean angular speed of revolution. If the undeter- 
mined scale factor be chosen so that Xi is the major semi-axis of the orbit of 
nil, then the mean angular velocity of revolution of the system is o>. 

The hypothesis is made that w* and m^, • • ., to„ are real positive numbers, 
and the problem is to find the number of real solutions of (6) for any value 
of n. For each of these solutions there is a six-fold infinity of collinear con- 
figurations, the six arbitrary parameters being the two Avhich define the plane 
of motion, the one which defines the orientation of the orbits in their plane, 
the one which determines the epochs at Avhich the bodies pass their apses, the 
one which determines the scale of the system, and finally the eccentricity 
of the orbits. 

3. Outline of the Method of Solution. The solution involves a 
mathematical induction and consists of the following steps : 
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Assumption (A) . To get the induction, it is assumed that for n = v the 
number of real solutions of (6) for x^, • • •, x^ is iV|„ whatever real positive 
values «* and nti, • • •, m, may have. It is known from the work of Lagrange 
that when i; := 3 we have i\^ — 3 = J 3 ! . 

Theorem (B). If to the system m^, • • •, to„ of positive masses an infi- 
nitesimal mass * m^^ i be added, then the whole number of real solutions is 

(v + 1) iv;. 

Theorem (0). As the infinitesimal mass m^^ increases continuously 
to any finite positive value whatever, the total number of real solutions re- 
mains precisely (i^ + 1) JV,. 

Conclusion (D) . By successive applications of theorems (B) and ( C) 
it is seen that the number of real solutions of (6) for n + i; + ju, is 

^r+^= ("+ M)(^ + M- 1) ••• (l'+ 2)(l/+ 1) iV",. 

When J/ = 3 it is known that JV3 = J3 !. Therefore 

^3 + ^ = J(/^ + 3)!. 
Let fi + 3 = n and we have 

(7) JSr,,=.inl 

To complete the demonstration of this conclusion it remains only to 
prove theorems (B) and (0). 

4. Proof of Theorem (B). When there are v finite bodies 
mj, . . ., m^ and the infinitesimal body m^+j, equations (6) become 



* The "infinitesimal body " Is In celestial mechanics a well established definite thing. In 
the first place it is absolutely zero. But obviously if this were all it would have no interest. 
It is used in association with certain differential equations which are derived in -writing 
them first for all the masses real and positive, then letting one of them (the infinitesimal mass) 
approach zero. The limit of this variable mass, zero, is the infinitesimal mass, and as this 
mass approaches zero, its motion satisfies the differential equations, which at the limit zero for 
the variable mass define the motion of the infinitesimal body. Since in this limiting process 
the differential equations remain determinate and in all respects regular the process is fully 
justified. It Is precisely in this sense that the infinitesimal mass is used in Theorem C-B) ; 
the equations used are the limits of equations for all masses finite as wi,- + 1 approaches 0. In 
Theorem (0) the massOT...|- 1 increases continuously and we follow the roots of the algebraic 
equations, at least as to their character. 
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The last column of these equations is zero because m^^i — 0, but is writ- 
ten for use in the proof of Theorem (O). Consequently the first v equations, 
which involve Xi, • • •, x^ alone as unknowns, are the equations defining the 
solutions when n = v. By {A) there are iVj, real solutions of these equations. 
Let any one of these solutions be x^ — x'J\ ■ ■ •, x^ = x'-'^. Then the last 
equation of (8) becomes 



^K+i = 



was. 



+1 



+ 



rrh(x^^^-o^ m^+i - 4''>) 



+ 



' 2r+l 



■ + 



+ 



m,(x^^i--a^ 



^Iv + l '2r+l 'rv+l 

The number of real solutions of this equation is required. 

Consider ^^+i as a function of x^^i. It is easily verified that 



(10) 



^^ + 1(4- co) = — 00 , 



limj, /™(^ 



e) = - 00, (y= 1, 

^^ + i(- 00) = + 00 . 






Since ^^+1 is finite and continuous except at x^^j = x^{\ • • •, x^^^ — x^°\ 
+ 00 , — 00 , it follows that there is an odd number of real solutions in each 
of the intervals — 00 to x^°\ where x^ is the smallest x'-f ; x'£? to x'-f, where x*^' 
and x^"' are any two x'j' which are adjacent ; and x<°' to + 00 , where x*°> is the 



largest x<j^. 
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But we find from (9) that 



which is negative except at a;^+i = Xj, . . ., x^^^ = a;^, where it is infinite. 
Therefore ^, + i is a decreasing monotonic function in each of the intervals, 
and consequently vanishes once, and but once, in each of them. Since there 
are k + 1 of these intervals there are, for each real solution of the first !< equa- 
tions of (8), precisely v ■{■ 1 real solutions of the last equation of (8). Since 
the first V equations have, by {A), N, real solutions, equations (8) altogether 
have precisely (v + 1) iVj, real solutions. This completes the demonstration 
of Theorem (J3) . 

4. Proof of Theorem (C). Let Xj = xf (J =1, • • -, v -\- l)be any 
one of the (v + 1)JV^ real solutions of equations (8) which are known to exist 
for m^+i = 0. It will be shown that as ot^+i increases continuously to any 
finite positive quantity whatever, the x'^p may be made to vary continuously so 
as always to satisfy equations (8) , and that during these variations the a;^°' remain 
distinct, finite, and real. From this it will follow that there are at least (v+ 1) JV^ 
real solutions of (8) for every set of finite positive values of wij, • • •, w^+j. It 
will also be shown that no new solutions can appear as m^^^ increases from 
zero to any finite value. Hence it will follow that the number of real solutions 
of (8) is exactly (v + 1) JV^ for all finite positive values of wii, • • •, m^^i. 

The roots of algebraic equations are continuous functions of the coeffi- 
cients so long as the roots are finite and the equations do not have indetermi- 
nate forms. Consequently the a;*"* are continuous functions of m^^i if no 
x''f becomes infinite and if no x^f> = x'-f. The real roots of algebraic equa- 
tions having real coefficients can disappear only by passing to infinity, or 
by an even number of real solutions becoming conjugate complex quantities 
in pairs. Therefore we have to determine (1) whether any finite x^f can 
become equal to any x^">, (2) whether any x'^' can become infinite, and 
(3) whether any two real solutions can become conjugate complex quan- 
tities for any finite positive values of mj, • • •, m^+i. 

(1). The masses m,, . . ., m^^i are by hypothesis all positive. Let 
the notation be chosen so that for any values of mi, • • •, to^+i for 
which the x'-'} are all distinct we have the inequalities x^ < x^l^ < • • • < x*; 
< ^y'\.\. Let us suppose that as some mass is changed we find x'-j^ — x\°^ 
approaching zero in such a way that x<j> and x^°> remain finite. That is, 
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ry, which occurs only in the expressions for </>j and <f>j, approaches zero. 
Suppose i < j. Then the term involving ry becomes negatively infinite 
in 4>i and positively infinite in <f>j. Consider <^,- = 0. Another »•,* must 
approach zero in order to restore the finite value of the function ^j, and 
the terra involving r^^. must become positively infinite as r,i approaches zero. 
Tlierefore we see from (8) that k < i. But ri^ enters besides only in </);t, 
and similar reasoning shows that r^-j, where I <h^ must also approach zero. 
In this manner we are driven to the conclusion finally that an rp^, where 
one of the subscripts is unity, approaches zero. Then consider </>i = 0. 
All its terms except — w* x^ are negative, and since one of its r^, tliat is 
Vpq, approaches zero, the first equation of (8) can not be satisfied. Con- 
sequently the original assumption that some r^ can approach zero for 
finite values of x^y\ • • • , aJ^ + i and finite positive values of mj, • • • , m^^j 
leads to an impossibility, and is therefore false. 

(2). Multiplying equations (8) by m^, m.i, • ■ •, ??i^_|_i respectively and 
adding, we get 

— afi(jn]Xy + riiiX^, + • • • + mjx^ + »"v+i ^r + i) — 0. 

It follows from this equation that no a;*f alone can become infinite, 
and that if one becomes positively infinite, then some other one must be- 
come negatively infinite. 

Suppose the notation is again chosen so that x*}' < a; 2' < • • • < x'°* < x^'^. 1. 
Then if any x'-f becomes negatively infinite x',' must do tlie same, and from 
the equation above it follows that arj,':|.i must become positively infinite. Now 
suppose this to occur and consider the equation ^1 = 0. In order that 
this equation may remain satisfied, x'|^ must also become negatively infinite 
in such a way that x^^ — x'-^ shall approach zero. But now it follows from 

<^2 = 0, since — (o^ x^°^ and ^^ — ^ are both po^tive, that x'j* must 

'"is 

also become negatively infinite in such a way that x*-"^ — x'-"^ shall approach zero. 
It follows similarly from ^3 = that x'p must become negatively infinite in 
such a way that x^"* — x^"* shall approach zero. This reasoning continues until 
it is found that x'f, • • •, xj,"^i must all become negatively infinite. But x'-'^i 
at least must become positively infinite, and we are thus led to a contra- 
diction. Similarly it can be proved that no x'-f can become positively 
infinite. 

In order to prove now that as wi^ + i approaches zero, equations (8) 
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remain determinate and their solutions finite, and that there are accordingly 
no others besides those obtained in Theorem (i?), consider a solution 
Xi, • • •, x,^i in which the Xj will be distinct for a set of positive values 
of mi, • • •, »i, + i, and then let m^ approach zero as a limit. 

In the first place, if Xf approaches neither Xi_i nor aj^^i as a finite limit 
as mi approaches zero as a limit, then by the reasoning of (1) and (2) above 
no Xj can approach any Xf, as a limit. 

In the second place, Xi can not approach ajj^.! as a finite limit as wij ap- 
proaches zero unless Xi_i also approaches Xj^j, for otherwise ^, = can 
not be satisfied. But if a;,_i approaches a;< + x as a limit as »ij approaches 
zero, then <j>i^i = and <f>i+i = can not be satisfied unless a;,_2 and Xj^j 
also approach Xf^i as a limit. This shifts the difficulty to <l>i-2 = and 
<}>i+i = 0, and so on until ^i = and ^„^i := are reached, which can 
not be satisfied under the hypothesis. 

In the third place, Xf can not become positively infinite as jra,- ap- 
proaches zero, for then ^i = can not be satisfied unless a;,_i becomes 
infinite in such a way that «,• — Xi_i approaches zero. Continuing through 
(^,_i = 0, etc., we are led to the conclusion that Xi, • • •, x^^i all become 
positively infinite ; but then the center of gravity equation can not be satis- 
fied. In a similar manner it can be proved that as wij approaches zero as a 
limit no Xj can approach any x^ or become infinite. Consequently the solu- 
tions all remain regular as mj approaches zero as a limit. 

(3). Since the solutions of (8) are continuous functions of m^^j, it fol- 
lows that no two solutions which are real for m^^i = can ever become con- 
jugate complex solutions for any real value of m^^i without having first 
become equal ; and similarly, no two solutions which are complex for 
m^^i = can ever become real for any real value of m^^j without 
having first become equal. Consequently if a multiple solution of (8) 
is impossible for every set of finite positive values of rrii, • ■ -, wi^ + i, 
it is impossible that any real solutions should disappear by becoming im- 
aginary, or that any imaginary solutions should become real. 

The conditions that x = x*"' shall be a multiple solution of f(x) = 
are that/ (x^"') = and /'(x*"') = 0. The corresponding conditions that a 
set of simultaneous algebraic equations shall have a multiple solution are 
that a set of values of the variables shall satisfy the equations and that 
the Jacobian of the functions with respect to the dependent variables shall 
vanish for the same set of values. That is, the conditions that Xj = xy\ 
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(y = 1, • • •, V + 1), shall be a multiple solution of (8) are that these 
values shall satisfy (8) and also the equation 



(11) 
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9^1 
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+ 1 



dxi dx2 



H. 



v+i 



dx, 



-+i 



= 0. 



Consider two solutions of a set of algebraic equations having real 
coeflScients. If they change from real to conjugate complex quantities or from 
conjugate complex to real quantities for a certain value of a continuously 
varying parameter, then for this value of the parameter, they are not only 
equal, but they are also real. Consequently we shall have occasion to exam- 
ine A only when all of its elements are real. We shall show that it cannot 
vanish for any finite set of real values of the Xj when m^, • • •, m^^^ 
are positive, and consequently that it can not vanish for any particular set 
which satisfies equations (8) . When this is established we shall have proved 
that all the solutions of (8) which are real for m^^., = remain real when 
m^ + i increases to any positive value, and those which are complex remain 
complex. 

From equations (8) and (11) we find 
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where 

~12 'iK+l 

2w,, 2^3 2m,,^i 

U 5 ... - 



'31 "23 '2i-|-l 



' + 1 



2 m, 2m, 



' + 11 '^4-1^ 



-0. 



If »«,-(- 1 = this determinant breaks up into the product of a determinant of 
the same type as (12) when all the m,- are distinct from zero and a factor es- 
sentially negative. Therefore it is sufficient to consider the genei^al case in 
which all the m^- are positive in examining whether A can vanish or not. 

Several properties of A are evident, (a) If the «th row be multiplied 
by mi (i — 1, • ■ -, V + 1), the determinant becomes symmetrical. (b) The 
sum of the elements in each row is — co^, from which it follows that the ex- 
pansion of the determinant contains w^ as a factor, (c) The expansion of the 
determinant contains (— 1)"+^ ft>2(''+i) as one of its terms, and since all the 
nil ^^^ positive and all the x, are real the sign of all the terms coming from 
the product of the elements of the main diagonal is (— 1)" + ^. 

When A is completely expanded all those terms not having the sign 
(— 1)"+^ are cancelled by terms coming from the product of the main 
diagonal elements, and since the term (— 1)"+' a)^(''+^) is certainly present 
the determinant can in no case be zero. The following demonstration of 
this fact was invented in 1907 by Mr. T. H. Hildebrandt, now of the Univer- 
sity of Michigan, as a class exercise.* 

Since the determinant contains w* as a factor every term in its expansion 
must depend upon at least one of the elements of the main diagonal. Fasten 
the attention upon any term of the expansion. It can be supposed without 
loss of generality that it depends upon the first main diagonal element. In 
the expansion of the determinant this element is multiplied by its minor ; 
consequently we must see if the minor can vanish. The minor is of the same 
form as the original determinant, and the sum of the elements of its ith row 

* An earlier proof was devised by the author, and still another jointly by Professor N. 
B. MacLeau, of the University of Manitoba, and Mr. E. J. Moultou, now of the University of 
Wisconsin. 
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is _ 0)2 _ 2mi/r|j. Consequently every term in the expansion of the minor 
will contain at least one of the — aP— 2«ii/rjj as a factor. But these elements 
appear only in the main diagonal of the minor. Hence all terms in the ex- 
pansion of the minor which do not vanish depend upon at least one element 
of the main diagonal. In considering our particular term it may be supposed 
without loss of generality to depend upon the first main diagonal element of 
the minor. In the expansion of the original determinant the product of these 
two diagonal elements will be multiplied by the co-factor of the minor of the 
second order of which they are the main diagonal. This co-factor has prop- 
erties similar to those of the first minor just considered, and in the same way 
we prove that at least one of its diagonal elements must be involved in the 
term in question. That is, the term under consideration depends upon at 
least three elements of the main diagonal. Continuing in this manner we 
prove that any term in the final expansion difierent from zero depends upon 
all the elements of the main diagonal which are all of the same sign in every 
one of their terms. Consequently all the terms which do not cancel out in 
the expansion of the determinant have the sign (— 1)" + '. And we have seen 
that there is at least one such term, viz., (— l)''+^ft)*('' + i). Therefore the deter- 
minant can not only never vanish, but it can never become less than fi)^<'' + ^) in 
numerical value. 

Since A can never vanish for real finite xf when all the m^ are real and 
positive, it follows that no real solutions can ever be lost or gained as the m^ 
vary, and therefore that the number of real solutions of (8) is always 
exactly {v + \) N^ = ^ {v + \) \ . 

5. Computation of the Solutions of Equations (6). There are 
well known methods of finding the roots of a single numerical algebraic equation 
of high degree, but they are not readily applicable to simultaneous equations of 
high degree. However, when the order of the masses has been chosen, equations 
(6) will become polynomials in x,, • • ., a;„ after they have been cleared of 
fractions. Then by rational processes n — 1 of the x^ can be eliminated from 
these equations giving a single equation in the remaining unknown. The 
solutions of this equation can be found by the usual methods, and the results 
can be used to eliminate one unknown. By repeated applications of this 
process to the successively reduced equations, the solutions can all be found. 
The single one satisfying the conditions of reality of Xj, • • ., x,, and their order 
relations is the one desired. 

The solutions can also be found by a method closely related to that by 
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means of which their existence was proved above. Suppose for w^ = mY^, 
(^ = 1, • • •, w), a solution a;,- = xj"'' of equations (6) is known. The wi|°^ are 
supposed to be zero or positive. Suppose it is desired to find the correspond- 
ing solution, that is, the one in which the masses are arranged on the line in 
the same order, for m^ = m^"> + ^f. Let the corresponding set of the Xi satis- 
fying (6) be Xf = x\'''> + fj, where the ^^ are functions of /^i, • • •, n„ to be de- 
termined. Substituting a;,- = a;J°> + f,- and m^ — wi/"* + /t,- in (6), making use 
of the notation of (8), expanding as power series in the fjand fii (which is 
always possible since we have shown no x\°'' can become infinite and no x\°'' can 
equal any xf^), and remembering that x,- = x^p is a solution of (6) for 
nii = m\°'', we have 

2 (>4>\ fc , v^ JL- "V ^^^ p\ 'V ^ 

i = 1 ■' ! — •> i — 1 J I 1 = 1 ■' 



-M. 



"p 



(13) 



2 £■«. ^ S 



./ = 1 



dxj 



2 



8x,- 



i, 



= -E 



i=x 



9^ 



f'P 



where 



2 



_9^ 

dxj 



h 



are the symbolic powers used in connection with 



the power series expansions of functions of several variables. 

The determinant of the terms of the first degree in the f^- in equations (13) 
is the A of equation (11) which we have proved can never vanish in this prob- 
lem. Therefore equations (13) can be uniquely solved for fj, • • •, f„ as power 
series in /itj, • • •, /t„ by the method of undetermined coefficients, and these 
series will converge for |;itj| >0 but sufficiently small. Suppose they con- 
verge if I ^j I S r. Keeping the /t,- within this limit a solution Xi = x\^'> is 
computed. Then this can be used as a starting point for a second application 
of the process, which can be repeated as many times as may be desired. 

Hence, to find the solution in which the bodies mj, • • •, m^ have any 
finite positive values and lie in a determined order on the line, we may start 
with mi, mj, and m^ and solve the Lagrangian quintic* which defines their dis- 

♦Tisserand's Mecanique Celeste, vol. 1, p. 155, or Moulton's Introduction to Celestial Me- 
chanics, p. 216. 
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tribution on the line. Then we add an infinitesimal body m^ and find its 
position by solving (9) in which now v = S. Then we increase the infinitesi- 
mal mass, step by step, to the required finite value, and compute the corres- 
ponding ajj, • • ., Xi. Making use of the fact that for |f,| < /a we have for all 
values of i 

\4>i\ <-4, ASa)«», 

where ^ is a number depending on p, co^, rui,- ■ ■, nin, it can be proved that by 
this process any finite value m^ can be reached in a finite number of steps. 
After the required value of nii has been reached the process can be repeated 
for wjg, etc., to any finite number of bodies. Notwithstanding the fact that 
this would be very laborious if the number of bodies was large, yet we must 
regard the problem as completely solved both theoretically and practically. 

II. Determination of the Masses when the Positions are Given. 

6. Determination of tlie Masses when n is Even. Suppose m^ 
and the n distinct points on a line, Xi, • • •, x„, are given, and consider the 
problem of determining mj, • • •, m„ so that the straight line solutions shall 
exist. There will be no loss of generality in selecting the notation so that 
Xi <x.i< ■ ■ -KX^. With this choice of notation equations (6), which are the 
necessary and sufficient conditions for the solutions, become 






+ 


m.i m^ 
rL 'is 


+ • • 


n n - 1 'in 


= - 0)2 Xy, 




+ 


+ ? 


-+ • • 


• ■>" „2 1 ,.2 

' 2 n — 1 '2n 


= — a>^ Xj, 



(14) 



m, m^ mg , »h, _ .^ 

. — — -J ... — u + — jj = —CO x„_i, 

n—ln 'n-li 'Ti— 1 3 '"»— !?» 

9W, m.2 _ ^ _ »»« 

' n\ ' «2 ' «3 ' nn 



+ = - ft)« a;„. 
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The m{ enter these equations linearly and are therefore uniquely 
determined if the determinant 



(15) 



Z) = 



'21 
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'12 
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*"«— 1 2 
- 1 



'»2 



K-IS 
- 1 



'jl3 



^«« + l 



' n—\ n 





is distinct from zero. This is a skew-symmetric determinant, and when n is 
even it is the square of an associated Pfaffian, and, therefore, is not in general 
zero. Therefore if n is even the masses are in general uniquely determined 
when 0)^ and Xj, • • •, x^ are given, though it should be noted that they are not 
necessarily all positive. When they are negative they have no physical in- 
terpretation. 

7. Determination of the Masses when n is odd. In this case 
the skew-symmetric determinant is identically zero, but its first minors of the 
main diagonal elements, being skew-symmetrical determinants of even order, 
are in general all distinct from zero ; consequently the x^ must satisfy one rela- 
tion in order that equations (14) shall be consistent. To get this relation add 
the equation 

niyXi + in.2X^ + . . . 4- m„ a;„ = 0, 



which is a consequence of (14), to the set of equations. In order that they 
shall be consistent their eliminant 
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must vanish. This is also a skew-symmetrical determinant and is the square 
of the PfaflBan 

F = \ Xi x^ • • • x„ 
J_ 1 

'19. /1„ 



'i2 



(16) 



= 0. 



Equation (16) can be found also bj solving any n — 1 equations of (14) 
for the corresponding TOj-and substituting the solutions in the remaining one. 
The result is a sum of determinants which can be shown to be the expansion 
of F multiplied by the square root of the determinant of the coefficients of 
the n — 1 OTj in the equations used. 

When i?'= is satisfied by «i, • • •, x„, equations (14) are consistent. 
Then, after any ot^ has been chosen arbitrarily, the corresponding ?i — 1 equa- 
tions can in general be solved uniquely for the remaining rrij, and the unused 
equation will be satisfied because F = 0. 

8. Discussion of Case n = 3. When n = 3 the determinant D 
becomes 

1 1 



and the Pfaffian F is 
(17) 



(^12 »-23 '•is)* ('-12 '"23 '•] 
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= 0, 



rT_ a^i 
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It will now be shown that when any two of x^, x^, x-^ are chosen so as to 
satisfy the conditions Xj < Xj < x.^, the third is uniquely determined by (17) 
and these inequalities. From the fact that in this case r^^ > r^^, »'i3 > r^ it 
follows that if Xj is positive then — x^/rls + x^/r'^ is positive, and therefore 
that Xi must be negative in order that (17) may be satisfied. If x^ is nega- 
tive Xj, being less, must also be negative. That is, x, is necessarily negative, 
and similarly x^ is necessarily positive. 

Suppose Xj and X3 are chosen and consider i^ as a function of Xj. Then 
we have at once 

,,„. lira i^(xi)=-oo. Mm F{x^-e) = + K, ^_ i 2x2 2X3 

j;, = — 00 0X1 > 23 "is ' 12 

dF 
From the inequalities Xg < X3 and r^ < rjg it follows that ^ — is positive 

cXj 

for Xi < Xj. Therefore there is one and but one solution of (17) for x, < Xj 
when Xj and a positive X3 are chosen. By symmetry there is but one solution 
of (17) for Xg > x-i when a negative Xx and x^ are chosen. 

When Xi and X3 are chosen respectively negative and positive but 
otherwise arbitrarily we have, considering i^ as a function of Xj, 

,,., limi^(Xi + e) = 4- 00 , lim i?'(x3 - €) = - 00 , ^'_!^i_ J__^^ 
(19) '=0 •=« 0x2 rl, r^, r^, < "' 

Therefore there is one and but one solution of (17) for X2 satisfying the ine- 
qualities Xi < Xj < X3. 

Suppose a negative Xj, a positive X3, and m^ are given arbitrarily and 
that x, is defined by (17). From (14) we have 



(20) 



Wli = - J-?3 



»l3 = - »*I3 



zr - «'a;3 

I-M2 



Since x^ is negative and X3 positive it follows from (20) that mj can be 
taken so small that mi and m^ are positive. If — Xj is not equal to X3, then 
a positive m^ can be determined so that nti and m^ shall both be positive, one 
positive and one negative, or both negative ; but wij, m.^, and m^ can not all be 
negative. 
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